Homologous mutations in human ß, embryonic, and perinatal muscle myosins have divergent effects on molecular power generation.

Mutations at a highly conserved homologous residue in three closely related muscle myosins cause three distinct diseases involving muscle defects: R671C in ß-cardiac myosin causes hypertrophic cardiomyopathy, R672C and R672H in embryonic skeletal myosin cause Freeman-Sheldon syndrome, and R674Q in perinatal skeletal myosin causes trismus-pseudocamptodactyly syndrome. It is not known whether their effects at the molecular level are similar to one another or correlate with disease phenotype and severity. To this end, we investigated the effects of the homologous mutations on key factors of molecular power production using recombinantly expressed human ß, embryonic, and perinatal myosin subfragment-1. We found large effects in the developmental myosins but minimal effects in ß myosin, and magnitude of changes correlated partially with clinical severity. The mutations in the developmental myosins dramatically decreased the step size and load-sensitive actin-detachment rate of single molecules measured by optical tweezers, in addition to decreasing overall enzymatic (ATPase) cycle rate. In contrast, the only measured effect of R671C in ß myosin was a larger step size. Our measurements of step size and bound times predicted velocities consistent with those measured in an in vitro motility assay. Finally, molecular dynamics simulations predicted that the arginine to cysteine mutation in embryonic, but not ß, myosin may reduce pre-powerstroke lever arm priming and ADP pocket opening, providing a possible structural mechanism consistent with the experimental observations. This paper presents direct comparisons of homologous mutations in several different myosin isoforms, whose divergent functional effects are a testament to myosin's highly allosteric nature.

Toward physics-based precision medicine: Exploiting protein dynamics to design new therapeutics and interpret variants.

The goal of precision medicine is to utilize our knowledge of the molecular causes of disease to better diagnose and treat patients. However, there is a substantial mismatch between the small number of food and drug administration (FDA)-approved drugs and annotated coding variants compared to the needs of precision medicine. This review introduces the concept of physics-based precision medicine, a scalable framework that promises to improve our understanding of sequence-function relationships and accelerate drug discovery. We show that accounting for the ensemble of structures a protein adopts in solution with computer simulations overcomes many of the limitations imposed by assuming a single protein structure. We highlight studies of protein dynamics and recent methods for the analysis of structural ensembles. These studies demonstrate that differences in conformational distributions predict functional differences within protein families and between variants. Thanks to new computational tools that are providing unprecedented access to protein structural ensembles, this insight may enable accurate predictions of variant pathogenicity for entire libraries of variants. We further show that explicitly accounting for protein ensembles, with methods like alchemical free energy calculations or docking to Markov state models, can uncover novel lead compounds. To conclude, we demonstrate that cryptic pockets, or cavities absent in experimental structures, provide an avenue to target proteins that are currently considered undruggable. Taken together, our review provides a roadmap for the field of protein science to accelerate precision medicine.

Homologous mutations in ß, embryonic, and perinatal muscle myosins have divergent effects on molecular power generation.

Mutations at a highly conserved homologous residue in three closely related muscle myosins cause three distinct diseases involving muscle defects: R671C in ß-cardiac myosin causes hypertrophic cardiomyopathy, R672C and R672H in embryonic skeletal myosin cause Freeman Sheldon syndrome, and R674Q in perinatal skeletal myosin causes trismus-pseudocamptodactyly syndrome. It is not known if their effects at the molecular level are similar to one another or correlate with disease phenotype and severity. To this end, we investigated the effects of the homologous mutations on key factors of molecular power production using recombinantly expressed human ß, embryonic, and perinatal myosin subfragment-1. We found large effects in the developmental myosins, with the most dramatic in perinatal, but minimal effects in ß myosin, and magnitude of changes correlated partially with clinical severity. The mutations in the developmental myosins dramatically decreased the step size and load-sensitive actin-detachment rate of single molecules measured by optical tweezers, in addition to decreasing ATPase cycle rate. In contrast, the only measured effect of R671C in ß myosin was a larger step size. Our measurements of step size and bound times predicted velocities consistent with those measured in an in vitro motility assay. Finally, molecular dynamics simulations predicted that the arginine to cysteine mutation in embryonic, but not ß, myosin may reduce pre-powerstroke lever arm priming and ADP pocket opening, providing a possible structural mechanism consistent with the experimental observations. This paper presents the first direct comparisons of homologous mutations in several different myosin isoforms, whose divergent functional effects are yet another testament to myosin's highly allosteric nature.

Divergent Molecular Phenotypes in Point Mutations at the Same Residue in Beta-Myosin Heavy Chain Lead to Distinct Cardiomyopathies.

In genetic cardiomyopathies, a frequently described phenomenon is how similar mutations in one protein can lead to discrete clinical phenotypes. One example is illustrated by two mutations in beta myosin heavy chain (ß-MHC) that are linked to hypertrophic cardiomyopathy (HCM) (Ile467Val, I467V) and left ventricular non-compaction (LVNC) (Ile467Thr, I467T). To investigate how these missense mutations lead to independent diseases, we studied the molecular effects of each mutation using recombinant human ß-MHC Subfragment 1 (S1) in in vitro assays. Both HCM-I467V and LVNC-I467T S1 mutations exhibited similar mechanochemical function, including unchanged ATPase and enhanced actin velocity but had opposing effects on the super-relaxed (SRX) state of myosin. HCM-I467V S1 showed a small reduction in the SRX state, shifting myosin to a more actin-available state that may lead to the "gain-of-function" phenotype commonly described in HCM. In contrast, LVNC-I467T significantly increased the population of myosin in the ultra-slow SRX state. Interestingly, molecular dynamics simulations reveal that I467T allosterically disrupts interactions between ADP and the nucleotide-binding pocket, which may result in an increased ADP release rate. This predicted change in ADP release rate may define the enhanced actin velocity measured in LVNC-I467T, but also describe the uncoupled mechanochemical function for this mutation where the enhanced ADP release rate may be sufficient to offset the increased SRX population of myosin. These contrasting molecular effects may lead to contractile dysregulation that initiates LVNC-associated signaling pathways that progress the phenotype. Together, analysis of these mutations provides evidence that phenotypic complexity originates at the molecular level and is critical to understanding disease progression and developing therapies.

Discovery of a cryptic pocket in the AI-predicted structure of PPM1D phosphatase explains the binding site and potency of its allosteric inhibitors.

Virtual screening is a widely used tool for drug discovery, but its predictive power can vary dramatically depending on how much structural data is available. In the best case, crystal structures of a ligand-bound protein can help find more potent ligands. However, virtual screens tend to be less predictive when only ligand-free crystal structures are available, and even less predictive if a homology model or other predicted structure must be used. Here, we explore the possibility that this situation can be improved by better accounting for protein dynamics, as simulations started from a single structure have a reasonable chance of sampling nearby structures that are more compatible with ligand binding. As a specific example, we consider the cancer drug target PPM1D/Wip1 phosphatase, a protein that lacks crystal structures. High-throughput screens have led to the discovery of several allosteric inhibitors of PPM1D, but their binding mode remains unknown. To enable further drug discovery efforts, we assessed the predictive power of an AlphaFold-predicted structure of PPM1D and a Markov state model (MSM) built from molecular dynamics simulations initiated from that structure. Our simulations reveal a cryptic pocket at the interface between two important structural elements, the flap and hinge regions. Using deep learning to predict the pose quality of each docked compound for the active site and cryptic pocket suggests that the inhibitors strongly prefer binding to the cryptic pocket, consistent with their allosteric effect. The predicted affinities for the dynamically uncovered cryptic pocket also recapitulate the relative potencies of the compounds (τb = 0.70) better than the predicted affinities for the static AlphaFold-predicted structure (τb = 0.42). Taken together, these results suggest that targeting the cryptic pocket is a good strategy for drugging PPM1D and, more generally, that conformations selected from simulation can improve virtual screening when limited structural data is available.

Accelerating Cryptic Pocket Discovery Using AlphaFold.

Cryptic pockets, or pockets absent in ligand-free, experimentally determined structures, hold great potential as drug targets. However, cryptic pocket openings are often beyond the reach of conventional biomolecular simulations because certain cryptic pocket openings involve slow motions. Here, we investigate whether AlphaFold can be used to accelerate cryptic pocket discovery either by generating structures with open pockets directly or generating structures with partially open pockets that can be used as starting points for simulations. We use AlphaFold to generate ensembles for 10 known cryptic pocket examples, including five that were deposited after AlphaFold's training data were extracted from the PDB. We find that in 6 out of 10 cases AlphaFold samples the open state. For plasmepsin II, an aspartic protease from the causative agent of malaria, AlphaFold only captures a partial pocket opening. As a result, we ran simulations from an ensemble of AlphaFold-generated structures and show that this strategy samples cryptic pocket opening, even though an equivalent amount of simulations launched from a ligand-free experimental structure fails to do so. Markov state models (MSMs) constructed from the AlphaFold-seeded simulations quickly yield a free energy landscape of cryptic pocket opening that is in good agreement with the same landscape generated with well-tempered metadynamics. Taken together, our results demonstrate that AlphaFold has a useful role to play in cryptic pocket discovery but that many cryptic pockets may remain difficult to sample using AlphaFold alone.

Discovery of a cryptic pocket in the AI-predicted structure of PPM1D phosphatase explains the binding site and potency of its allosteric inhibitors.

Virtual screening is a widely used tool for drug discovery, but its predictive power can vary dramatically depending on how much structural data is available. In the best case, crystal structures of a ligand-bound protein can help find more potent ligands. However, virtual screens tend to be less predictive when only ligand-free crystal structures are available, and even less predictive if a homology model or other predicted structure must be used. Here, we explore the possibility that this situation can be improved by better accounting for protein dynamics, as simulations started from a single structure have a reasonable chance of sampling nearby structures that are more compatible with ligand binding. As a specific example, we consider the cancer drug target PPM1D/Wip1 phosphatase, a protein that lacks crystal structures. High-throughput screens have led to the discovery of several allosteric inhibitors of PPM1D, but their binding mode remains unknown. To enable further drug discovery efforts, we assessed the predictive power of an AlphaFold-predicted structure of PPM1D and a Markov state model (MSM) built from molecular dynamics simulations initiated from that structure. Our simulations reveal a cryptic pocket at the interface between two important structural elements, the flap and hinge regions. Using deep learning to predict the pose quality of each docked compound for the active site and cryptic pocket suggests that the inhibitors strongly prefer binding to the cryptic pocket, consistent with their allosteric effect. The predicted affinities for the dynamically uncovered cryptic pocket also recapitulate the relative potencies of the compounds (τ b =0.70) better than the predicted affinities for the static AlphaFold-predicted structure (τ b =0.42). Taken together, these results suggest that targeting the cryptic pocket is a good strategy for drugging PPM1D and, more generally, that conformations selected from simulation can improve virtual screening when limited structural data is available.

Predicting locations of cryptic pockets from single protein structures using the PocketMiner graph neural network.

Cryptic pockets expand the scope of drug discovery by enabling targeting of proteins currently considered undruggable because they lack pockets in their ground state structures. However, identifying cryptic pockets is labor-intensive and slow. The ability to accurately and rapidly predict if and where cryptic pockets are likely to form from a structure would greatly accelerate the search for druggable pockets. Here, we present PocketMiner, a graph neural network trained to predict where pockets are likely to open in molecular dynamics simulations. Applying PocketMiner to single structures from a newly curated dataset of 39 experimentally confirmed cryptic pockets demonstrates that it accurately identifies cryptic pockets (ROC-AUC: 0.87) >1,000-fold faster than existing methods. We apply PocketMiner across the human proteome and show that predicted pockets open in simulations, suggesting that over half of proteins thought to lack pockets based on available structures likely contain cryptic pockets, vastly expanding the potentially druggable proteome.

Drug specificity and affinity are encoded in the probability of cryptic pocket opening in myosin motor domains.

The design of compounds that can discriminate between closely related target proteins remains a central challenge in drug discovery. Specific therapeutics targeting the highly conserved myosin motor family are urgently needed as mutations in at least six of its members cause numerous diseases. Allosteric modulators, like the myosin-II inhibitor blebbistatin, are a promising means to achieve specificity. However, it remains unclear why blebbistatin inhibits myosin-II motors with different potencies given that it binds at a highly conserved pocket that is always closed in blebbistatin-free experimental structures. We hypothesized that the probability of pocket opening is an important determinant of the potency of compounds like blebbistatin. To test this hypothesis, we used Markov state models (MSMs) built from over 2 ms of aggregate molecular dynamics simulations with explicit solvent. We find that blebbistatin's binding pocket readily opens in simulations of blebbistatin-sensitive myosin isoforms. Comparing these conformational ensembles reveals that the probability of pocket opening correctly identifies which isoforms are most sensitive to blebbistatin inhibition and that docking against MSMs quantitatively predicts blebbistatin binding affinities (R2=0.82). In a blind prediction for an isoform (Myh7b) whose blebbistatin sensitivity was unknown, we find good agreement between predicted and measured IC50s (0.67 µM vs. 0.36 µM). Therefore, we expect this framework to be useful for the development of novel specific drugs across numerous protein targets.

Functional divergence of the sarcomeric myosin, MYH7b, supports species-specific biological roles.

Myosin heavy chain 7b (MYH7b) is an evolutionarily ancient member of the sarcomeric myosin family, which typically supports striated muscle function. However, in mammals, alternative splicing prevents MYH7b protein production in cardiac and most skeletal muscles and limits expression to a subset of specialized muscles and certain nonmuscle environments. In contrast, MYH7b protein is abundant in python cardiac and skeletal muscles. Although the MYH7b expression pattern diverges in mammals versus reptiles, MYH7b shares high sequence identity across species. So, it remains unclear how mammalian MYH7b function may differ from that of other sarcomeric myosins and whether human and python MYH7b motor functions diverge as their expression patterns suggest. Thus, we generated recombinant human and python MYH7b protein and measured their motor properties to investigate any species-specific differences in activity. Our results reveal that despite having similar working strokes, the MYH7b isoforms have slower actin-activated ATPase cycles and actin sliding velocities than human cardiac ß-MyHC. Furthermore, python MYH7b is tuned to have slower motor activity than human MYH7b because of slower kinetics of the chemomechanical cycle. We found that the MYH7b isoforms adopt a higher proportion of myosin heads in the ultraslow, super-relaxed state compared with human cardiac ß-MyHC. These findings are supported by molecular dynamics simulations that predict MYH7b preferentially occupies myosin active site conformations similar to those observed in the structurally inactive state. Together, these results suggest that MYH7b is specialized for slow and energy-conserving motor activity and that differential tuning of MYH7b orthologs contributes to species-specific biological roles.

New Copper Alloys Used to Make Products Intended for Contact with Drinking Water.

This article presents the results of tests conducted as part of a research project with the primary objective of developing new copper alloys with limited lead content. The new group of materials were created in a production plant. As part of tests, a group of 22 alloys were selected for testing in castability, structural characteristics and hardness. Based on the test results obtained, the group of alloys under study was narrowed down to nine. The mechanical properties of these alloys were determined in static tensile tests as well as in uniaxial upsetting tests at elevated temperature, on the basis of which the group of alloys under investigation was further narrowed to three. Further studies involved technological verification of the application of these alloys under industrial conditions. These alloys were subject to numerical forging analyses, along with forging tests, under semi-industrial conditions, where the degree of filling of a die impression at a specific temperature was measured using an optic scanner. The quality of production of the obtained forgings was evaluated macroscopically with simultaneous observations of the microstructure.

SARS-CoV-2 simulations go exascale to predict dramatic spike opening and cryptic pockets across the proteome.

SARS-CoV-2 has intricate mechanisms for initiating infection, immune evasion/suppression and replication that depend on the structure and dynamics of its constituent proteins. Many protein structures have been solved, but far less is known about their relevant conformational changes. To address this challenge, over a million citizen scientists banded together through the Folding@home distributed computing project to create the first exascale computer and simulate 0.1 seconds of the viral proteome. Our adaptive sampling simulations predict dramatic opening of the apo spike complex, far beyond that seen experimentally, explaining and predicting the existence of 'cryptic' epitopes. Different spike variants modulate the probabilities of open versus closed structures, balancing receptor binding and immune evasion. We also discover dramatic conformational changes across the proteome, which reveal over 50 'cryptic' pockets that expand targeting options for the design of antivirals. All data and models are freely available online, providing a quantitative structural atlas.

Deep learning the structural determinants of protein biochemical properties by comparing structural ensembles with DiffNets.

Understanding the structural determinants of a protein's biochemical properties, such as activity and stability, is a major challenge in biology and medicine. Comparing computer simulations of protein variants with different biochemical properties is an increasingly powerful means to drive progress. However, success often hinges on dimensionality reduction algorithms for simplifying the complex ensemble of structures each variant adopts. Unfortunately, common algorithms rely on potentially misleading assumptions about what structural features are important, such as emphasizing larger geometric changes over smaller ones. Here we present DiffNets, self-supervised autoencoders that avoid such assumptions, and automatically identify the relevant features, by requiring that the low-dimensional representations they learn are sufficient to predict the biochemical differences between protein variants. For example, DiffNets automatically identify subtle structural signatures that predict the relative stabilities of ß-lactamase variants and duty ratios of myosin isoforms. DiffNets should also be applicable to understanding other perturbations, such as ligand binding.

SARS-CoV-2 Simulations Go Exascale to Capture Spike Opening and Reveal Cryptic Pockets Across the Proteome.

SARS-CoV-2 has intricate mechanisms for initiating infection, immune evasion/suppression, and replication, which depend on the structure and dynamics of its constituent proteins. Many protein structures have been solved, but far less is known about their relevant conformational changes. To address this challenge, over a million citizen scientists banded together through the Folding@home distributed computing project to create the first exascale computer and simulate an unprecedented 0.1 seconds of the viral proteome. Our simulations capture dramatic opening of the apo Spike complex, far beyond that seen experimentally, which explains and successfully predicts the existence of 'cryptic' epitopes. Different Spike homologues modulate the probabilities of open versus closed structures, balancing receptor binding and immune evasion. We also observe dramatic conformational changes across the proteome, which reveal over 50 'cryptic' pockets that expand targeting options for the design of antivirals. All data and models are freely available online, providing a quantitative structural atlas.

Conformational distributions of isolated myosin motor domains encode their mechanochemical properties.

Myosin motor domains perform an extraordinary diversity of biological functions despite sharing a common mechanochemical cycle. Motors are adapted to their function, in part, by tuning the thermodynamics and kinetics of steps in this cycle. However, it remains unclear how sequence encodes these differences, since biochemically distinct motors often have nearly indistinguishable crystal structures. We hypothesized that sequences produce distinct biochemical phenotypes by modulating the relative probabilities of an ensemble of conformations primed for different functional roles. To test this hypothesis, we modeled the distribution of conformations for 12 myosin motor domains by building Markov state models (MSMs) from an unprecedented two milliseconds of all-atom, explicit-solvent molecular dynamics simulations. Comparing motors reveals shifts in the balance between nucleotide-favorable and nucleotide-unfavorable P-loop conformations that predict experimentally measured duty ratios and ADP release rates better than sequence or individual structures. This result demonstrates the power of an ensemble perspective for interrogating sequence-function relationships.

Left hemisphere structural connectivity abnormality in pediatric hydrocephalus patients following surgery.

Neuroimaging research in surgically treated pediatric hydrocephalus patients remains challenging due to the artifact caused by programmable shunt. Our previous study has demonstrated significant alterations in the whole brain white matter structural connectivity based on diffusion tensor imaging (DTI) and graph theoretical analysis in children with hydrocephalus prior to surgery or in surgically treated children without programmable shunts. This study seeks to investigate the impact of brain injury on the topological features in the left hemisphere, contratelateral to the shunt placement, which will avoid the influence of shunt artifacts and makes further group comparisons feasible for children with programmable shunt valves. Three groups of children (34 in the control group, 12 in the 3-month post-surgery group, and 24 in the 12-month post-surgery group, age between 1 and 18 years) were included in the study. The structural connectivity data processing and analysis were performed based on DTI and graph theoretical analysis. Specific procedures were revised to include only left brain imaging data in normalization, parcellation, and fiber counting from DTI tractography. Our results showed that, when compared to controls, children with hydrocephalus in both the 3-month and 12-month post-surgery groups had significantly lower normalized clustering coefficient, lower small-worldness, and higher global efficiency (all p < 0.05, corrected). At a regional level, both patient groups showed significant alteration in one or more regional connectivity measures in a series of brain regions in the left hemisphere (8 and 10 regions in the 3-month post-surgery and the 12-month post-surgery group, respectively, all p < 0.05, corrected). No significant correlation was found between any of the global or regional measures and the contemporaneous neuropsychological outcomes [the General Adaptive Composite (GAC) from the Adaptive Behavior Assessment System, Second Edition (ABAS-II)]. However, one global network measure (global efficiency) and two regional network measures in the insula (local efficiency and between centrality) tested at 3-month post-surgery were found to correlate with GAC score tested at 12-month post-surgery with statistical significance (all p < 0.05, corrected). Our data showed that the structural connectivity analysis based on DTI and graph theory was sensitive in detecting both global and regional network abnormality when the analysis was conducted in the left hemisphere only. This approach provides a new avenue enabling the application of advanced neuroimaging analysis methods in quantifying brain damage in children with hydrocephalus surgically treated with programmable shunts.